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Abstract 

The vacuum diagram is calculated at second order for theories with self-interacting massless 
fields in the framework of finite causal perturbation theory. It is pointed out that the infrared 
behaviour of the vacuum diagram leads to unstable Fock vacua for QCD or massless QED, but not 
for quantum gravity. Therefore a radical rearrangement of the physical system must take place for 
such theories. Conversely, stability of the Fock vacuum for massless interacting fields is another 
hint at the possibility that quantum gravity should be treated as an effective theory. 
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1 Introduction 



Soon after the beginnings of Quantum Field Theory (QFT), it turned out that the perturbative 
expansion of several quantities like the electron self-energy are plagued by divergences of different 
kind. These divergences in QFT can be classified into four types, namely infrared, ultraviolet (UV), 
infinite volume and particle number divergences. A large amount of work has been done to overcome 
these problems, but even after the construction of a renomalized perturbation expansions based on 
the Lagrangean formalism for interacting fields or the scattering matrix (S-matrix), little is known 
until today about the convergence properties of the perturbation expansion. 

The first attempt to go beyond the widely used Lagrangean framework can be traced back to W. 
Heisenberg in 1943 jj], where he came to the conclusion that the basic observable is the S'-matrix, 
which describes actually the quantities that can be measured in an experiment. Heisenberg even 
proposed to construct a theory directly in terms of the elements of the S'-matrix, eliminating the 
notion of fields and the hypothesis of adiabatic switching of the interaction. However, the complete 
banishing of local quantities from the theory turned out to be too radical. 

A very natural and rigorous approach to axiomatic perturbation theory which has been applied 
successfully to all relevant interactions of the Standard model |2|, ||, |[] was provided by H. Epstein and 
V. Glaser in 1971 f|, |j 0]. Their method, often called finite causal perturbation theory or FCPT for 
short in this paper, avoids UV-divergences from the start by defining mathematically correct time- 
ordered products for distributions. The resulting distributions are smeared out by test functions, 
avoiding provisionally volume and infrared divergences at once. In FCPT, the ansatz for the S-matrix 
as a power series in the coupling constant is crucial, namely S is considered as a sum of smeared 
operator-valued distributions 

00 1 f 

S(g) = 1 + ^2 — /d 4 zi • • .d 4 x n T n (xi, . . . ,x n )g(xi) ■ . . . ■ g(x n ) , (1) 

n=l 

where the Schwartz test function g E 6*(R 4 ) switches the interaction and provides the infrared cutoff. 
The basic formulation of causality in FCPT, which had already been used by Bogoliubov et al. Q, is 

S(gi + g2) = S(g 2 )S( gi ) , (2) 

if the support of g2(x) is later than suppgi in some Lorentz frame. We denote this by supp gi > suppgi. 
The condition allows the construction of the n-point distributions T n as a well-defined 'renormalized' 
time-ordered product expressed in terms of Wick monomials of free fields :0(x\, . . . ,x n ): 

T n {x\, . . . , Xn) — ^ ' • 0{x\ , . . . , Xn) . t n (x\ X n , . . . , X n — 1 X n ) ■ (3) 
O 

Possible ambiguities in the decomposition of T n according to (||) will play no role for us in the sequel. 
The t n are C-number distributions. T n is constructed inductively from the first order T\(x), which 
describes the interaction among the quantum fields, and from the lower orders Tj, j = 2, . . . , n — 1 by 
means of Poincare covariance and causality. The inductive construction of the n-point distributions 
T n can be considered as the main technical drawback of the theory, since all lower orders Xi, ...T n _i 
must be calculated first in order to construct T n . But a naive definition of T n like 

£ e« - 4 2 ) • ... • 0(4^ - xIjt^) ■ ... • t^j (4) 

TT 
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where the sum runs over all n! permutations, is not well-defined, since it contains the product of 
Heaviside distributions with other singular distributions. This error leads to the well-known UV 
divergences in the calculation of Feynman diagrams. We illustrate this fact by the following very 
simple example in one dimension: It is of course meaningless to multiply a Heaviside distribution 0(i) 
with a 5-distribution 5(t) or even its derivative 5'(t). The Fourier transforms are 

©M = -4=^ ■ H») = J %= ■ (5) 
V27T lo -iO V2vr 

The ill-defined product (Q5')(t) goes over into a non-existing convolution by Fourier transform 

1 i r+oo , , , 

m > m ^- ( e , = 1 ^f_^^— (6 ) 

which is definitely 'UV-divergent'. 

T2(x, y) contains typically tree (one contraction), loop (two contractions) and vacuum graph (three 
contractions) contributions. Due to the presence of normal ordering, tadpole diagrams do not show 
up. In this paper, we will concentrate on the vacuum graph, which has been considered already 
in H in a preliminary way. Since this graph contains no external fields, we can directly identify 
T2 ac (xi,X2) = t2 ac (xi — X2). Usually, vacuum bubbles are neglected by simply absorbing them into 
a phase of the S'-matrix. In FCPT, this is not necessary, since their contribution to the S-matrix 
vanishes in the adiabatic limit for theories like ordinary QED with a massive electron. But this is not 
always the case for theories with self-interacting massless fields like QCD. 

2 Preliminary remarks 

The first order interaction T\(x) ~ i : Ci n t : is usually motivated from a classical Lagrangean interac- 
tion density. For QED, T\(x) is given by means of the free electron and photon field 

Ti(x) = ie : ~$(x)^y(x) : A^x) . (7) 

For the purpose of this paper, we outline briefly the main steps in the construction of T2(x, y) from a 
given first order interaction. In FCPT, one constructs first the operator valued distributions 

R' 2 (x,y) = -T x (y)ri(x) , (8) 
A' 2 (x,y) = -T^T^y) , and (9) 
D 2 (x,y) = R! 2 (x,y)-A 2 (x,y) = [T l {x),T 1 {y)} . (10) 

A meaningful theory ensures causality conditions (in a distributive sense) like (@), which implies also 

[T 1 (x),T 1 (y)}=0 for (x - yf < (11) 

such that 

sup V D2(x,y)Q{(x,y)e-R 8 \(x-y) 2 >0} = V2 . (12) 

This property of D2 must be checked first before one may start the inductive construction of the 
S'-matrix. From D2(x,y) we construct the Lorentz covariant retarded distribution R 2 (x,y), which 
coincides with D 2 (x,y) in the sense that (R 2 ,f) = (D 2 ,(p) for all test functions <p £ 5(R 8 ) with 
(compact) suppy? C {(x,y) G R 8 |(x° — y°) > 0} and (R 2 ,ip) = for all ip with supp <p C {(x,y) G 
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R 8 |(x° — y°) < 0}. Such a distribution can always be found, but it is only denned up to local 
operator-valued distributions with support on the complete diagonal A2 = {(x,y) £ R 8 |x = y}. 
Symmetry and renormalizability considerations further restrict these ambiguous terms. Obviously, 
T2(x,y) = R 2 (x,y) — R' 2 (x,y) is then a time-ordered product, since by construction 

M*,y)- j Tl{y)Tl[x) . x o_ y o <0 ■ 

We pursue the strategy described above for the calculation of the vacuum diagram described by the 
C-number distribution t 2 ac (x — y). 



3 Calculation of the vacuum diagram 

We calculate the vacuum diagram for the example of pure QCD, where the first order self-interaction 
is described by the help of 5f7(3) structure constants f abc 

Tx(x) = igf ahc :A^{x)A vh (x)d"A»(x)-. 

-igfabc ■ A^(x)u b (x)d^u c (x) : . (14) 

Here, A^ a are the free gauge potentials, satisfying the commutation relations 

[A^(x),A^ ± \y)]=iS ab g^Dt(x-y) , (15) 

where A^ are the emission and absorption parts of A and the zero-mass Pauli-Jordan distribu- 
tions defined by 

D±{x) = J d 4 P e(±p°)5(p 2 )e-^ . (16) 

The free ghosts are Fermi fields and satisfy the anticommutation relations 

{ U ^( x ),u { ±\y)} = -i5 ab D±(x-y) , (17) 

and all other anticommutators vanish. The quartic gluon interaction is missing in T\ since it is enforced 
by gauge invariance as a local counterterm at second order in FCPT Q. For the sake of convenience, 
we will use g for the QCD coupling constant as well as for test functions, since a mix up is impossible. 
As described in the previous section, we calculate first the vacuum expectation value 

B! 2 vac (x,y) = r'r c (x -y) = -< 0|r 1 (y)T 1 (x)|^ > . (18) 

The calculation is straightforward, therefore we present the result which can be decomposed in a 
purely gluonic part and a ghost part: 

r' 2 vac ' 9 \x -y) = -3ig 2 N(N 2 - l)D^{x - y)cf u D^x - y)d^(x - y) , (19) 

r' 2 vac ' gh (x -y) = +ig 2 N(N 2 - l)D (x - y)d x u D^(x - y)d^(x - y) , (20) 

where we have used that the fully antisymmetric structure constants satisfy fabcfabc = N(N 2 — 1) for 
SU{N). 
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It is difficult to give a meaning to the product of three singular objects in x-space, but r'^" 10 is of 
course well-defined via the impulse space, where the products in x-space go over into convolutions. 
The C-number distribution we consider now is therefore (omitting temporarily some prefactors) 

r'r%x) = -£>o (x)a M £?o (x)^D^(x) = -\d^{D^{x)f , (21) 
frW^ft-^o^^) . (22) 

The preliminary result 

(D * Do )(p) = -^e(p 2 M-p°) (23) 
can be readily obtained. Then we have 

24(2vr) 

If p is spacelike, i.e. p 2 < 0, then we can choose a special Lorentz frame where p = (0,p), and it follows 
that r£ ac (p) = due to the term @(q°)0(-q°) in (|J). For timelike p we choose a Lorentz frame 
where p = (p°,0). The four-dimensional integral can be restricted to three dimensions because of the 
^-distribution in (p — q) 2 . Taking the Heaviside distributions in ( g4j) into account, 

(P " q) 2 =vl~ 2p°q° + ql -q 2 = (25) 
implies q° = p° + | q\ and the integral goes over into 

ip\ f d 3 q 



rT c (v) = ^S-^b / Ai((p-«) 2 )e(?°-P°)e( ? 2 )e(- ? ) . (24) 



^° C (P°>0) = ^e/||0(Po + 2p O k1)e(-|g1-p ) 



,2 



*fo e(-p°)Po • ( 26 ) 
192(2vr) 5 v f ^° v ; 



The final result is by Lorentz covariance 

fr(p)= 192p e( " P ° )e(p2)(|,2)2 ' (27) 
From symmetry considerations we have automatically 

aT C {p) = rTi-p) (28) 

and therefore 

d v 2 ac (p) = -^^sgn^eb 2 )^ 2 ) 2 . (29) 

The calculation of the retarded part r™ c of d™ is demonstrated for an analogous case in [jl(]] in full 
detail. The result turns out to be 
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and the final result for the vacuum diagram is 

h [P) - 96(2vr)6 {P } l0g M 2 " 

Scale invariance of the causal distribution is destroyed by the time-ordering procedure. The parameter 
M 2 remains undefined. Changing the value of M 2 changes t v 2 ac only by local terms in x-space, which 
are not relevant for the following discussion. 

4 Volume divergence 

A problem arises when we calculate the vacuum transition amplitude at second order 

<n\S 2 (n)\n> = ^Jd A xd 4 yg(x)t v 2 ac (x-y)g(y) 

= i(2vr) 2 Jd 4 pg(-p)ir c (p)g(p) (32) 

in the adiabatic limit g — ► 1. Since 1 5(R 4 ), we have to specify the limiting process. Usually one 
chooses a fixed test function go G 5(R 4 ) with go(0) = 1 and performs a scaling limit, i.e. for e — > 
we have g e (x) = go(ex) — > 1. This implies 

<n\S 2 (n)\n> = ^Jd 4 xd 4 yg (ex)t v 2 ac (x-y)g (ey) 

= 2? / d 4 *A <7oW c ((z " y)/ £ )^(y) 

= di P9o(-P)i V 2 ac (ep)g (p) • (33) 



This expression diverges in the adiabatic limit, since 

M 2 



hi 



d 4 P go(-p)(ep) 4 log {£2p2 ^ l0) So( P ) - (34) 



J d 4 pg (-p)p 4 log y'l oM + log(£ 2 ) / A5o(p)p 4 5o(-p) (35) 
due to the logarithmic term in e. 

5 Conclusion 

The Fock vacuum is unstable against the interaction of the massless gluon fields. This is a consequence 
of the factor p 4 in the vacuum diagram, which expresses simply its power-counting degree u = 4. The 
instability has nothing to do with the non-abelian character of QCD, since the same phenomenon also 
occurs in massless QED. There, t 2 ac has exactly the same form as in QCD up to a negative prefactor 

^log^feLM • (36) 



1 24(2vr) 6 ^ b M' 2 

It is tempting to look for a cancellation of the QCD vacuum divergence through massless fermions, 
but in view of the present properties of QCD (confinement), this is probably not even desirable. 
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It is interesting to observe that in quantum gravity, the situation is remedied by the non-renormalizable 
coupling that leads to a factor p e in ti which cancels the logarithmic divergence. The same is true for 
the Euler effective action for photon-photon scattering 



Hn 2 - - -, 
Ti(x) = ^[(E 2 -B 2 f + 7(EBf] . (37) 



One might therefore argue for QCD that the volume divergence encountered in eq. fl3q ) must be 
compensated by the interaction, which signals a radical rearrangement of the system as in statistical 
mechanics. Such a rearrangement is well-known from (l+l)-dimensional QFT models, especially for 
the Schwinger model [|ll], 12, 13 1. 



Our result shows also that it is really a fair assumption that quantum gravity should be treated as 
an effective theory |I3j|, which incorporates the 'true' physical vacuum. The price one has to pay for 
perturbative vacuum stability is non-renormalizability, since we do not know the underlying theory of 
quantum gravity 
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